Field theoretic calculation of scalar turbulence 
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The cascade rate of passive scalar and Bachelor's constant in scalar turbulence are calculated 
using the flux formula. This calculation is done to first order in perturbation series. Batchelor's 
constant in three dimension is found to be approximately 1.25. In higher dimension, the constant 
increases as d 1//3 . 
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Perturbative field-theoretic techniques have been very useful in turbulence research. One of the celebrated field- 
theoretic method, renormalization groups (RG), has been applied to fluid turbulence [1, 2, 3], scalar turbulence [2, 4], 
MHD turbulence [5] etc. In RG analysis, one can calculate the renormalized parameters at large length scales. In 
addition to RG, one can also apply the field-theoretic techniques to calculate turbulent cascade rates [5, 6]. In this 
paper we will calculate the cascade rates of passive scalar using perturbative technique. From this calculation we can 
also calculate Batchelor's constant, which is very important for large-eddy simulations. 

The study of passive scalar is one of the important areas in turbulence research. It finds application in evolution 
of temperature field, pollution diffusion, etc. The phenomenology of passive scalar is well developed [6], and their 
predictions arc in agreement with the experimental results. According to the phenomenology, the energy spectrum 
of both velocity field u and scalar field tp in the inertial-convective range are proportional to k~ 5 ^ 3 . Note that in the 
inertial-convective range both the nonlinear terms u • Vu and u • V"0 dominate the viscous term. However, there exist 
Ch ' two other ranges depending on the value of Prandtl number (the ratio of viscosity and diffusivity) . In this paper we 
| will only focus on inertial-convective range. 

i Regarding the calculation of renormalized viscosity and diffusivity for passive scalar admixture, Yakhot and Orszag 
J> ' [2] adopted e-expansion, while Zhou and Vahala [4], and Lin et al.'s [7] procedure is recursive based on the original 
7— I , idea of McComb and his group ([3] and reference therein). Adzhemyan et al. [8] used De Dominicis and Martin's [9] 
procedure for fluid turbulence to passive scalars and computed the renormalized parameters. Earlier, Wyld [10] had 
given a perturbative expansion of Navier-Stokes equation. Canuto and Dubovikov [11], and Canuto et al. [12] started 
with Wyld's formalism and computed the renormalized diffusivity for passive scalar; they also computed Batchelor's 
constant. 

Turbulence cascade rates or fluxes play important role in turbulence calculations. It is a measure of transfer 
of a certain quantity from inside of a wavenumber sphere to the outside wavenumber sphere. In fluid turbulence, 
Kraichnan [13] applied direct interaction approximation and calculated the flux. Later, the cascade rates have been 
calculated by many researchers using various techniques, e.g. Eddy-Damped-Quasi-Normal-Markovian (EDQNM) 
closure scheme, RG, etc. Here we are interested in the cascade rate of passive scalar. This cascade rate quantifies 
how scalar fluctuations at large length-scales diffuse to small length-scales. 

In this paper we apply perturbative techniques to calculate the cascade rate of passive scalar. In our scheme the 
cascade rate of passive scalar is calculated using the flux formula and the renormalized parameters. In section 2 we 
recapitulate the earlier RG calculation [4] and extend their results to higher dimensions. In the subsequent section 
we apply the perturbative technique and calculate the cascade rate in the inertial-convective range to first order. The 
final expression involves energy spectrum for which we substitute /c~ 5 / 3 obtained from the phenomenology. From 
this procedure we also calculate Batchelor's constant. We have extended our calculations to higher space dimensions, 
because higher-dimensional field theory usually provide important insights into the nature of nonlinear interactions 

The outline of the paper is as follows: in section 2 we provide the definitions and recapitulation of the renormalization 
procedure for passive scalar. In section 3, we carry out the calculation of flux of passive scalar and Batchelor's constant. 
Section 4 contains conclusions. 



II. RENORMALIZATION OF VISCOSITY AND DIFFUSIVITY REVISITED 



Earlier calculations of renormalization in scalar turbulence have been carried out by Yakhot and Orszag [2] , Zhou 
and Vahala [4], and Lin et al. [7]. In this section we recapitulate very briefly Zhou and Vahala's calculation for passive 
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scalar and extend their results to higher dimensions. Zhou and Vahala's calculation is based on recursive scheme 
proposed by McComb and his coworker ([3] and references therein). The equations for the velocity u and passive 
scalar ip fields in Fourier space are 



{—iu + vk 2 ) Ui(k) = -^fijm(k) J dpuj(p)u m (k - p) (1) 

{ -i. + ^)m - -ik.jdpu^ik-p) (2) 

PijmQ^) = kjPi m (\t) — k m Pij(\i); (3) 

PimQ*-) — $im j (4) 

k = (k,w); (5) 
dp = dpdcj/(2Tr) d+1 . (6) 

Here v and k are the viscosity and diffusivity respectively, p is the fluid pressure, and d is the space dimension. We 
have assumed that the flow is incompressible, i.e. kiUi(k) = 0. 

In the recursive RG procedure the wavenumber range (fcjv, k ) is divided logarithmically into N shells. The effective 
parameters are obtained by eliminating the high wavenumber shells iteratively. We denote the higher wavenumber 
shells by k y and the remaining wavenumber region by fc < . In this procedure the field variables uf(k) and ^(k) are 
assumed to be gaussian with zero mean, and 

(u>(p)u>(q)) = P l3 (p)C u (p)S(P + q) (7) 
(V^ = C^(p)S(p + q) (8) 

where C u {p) and C^(p) are velocity and scalar correlation functions respectively. 

If we denote v^ n ) and «(„) as viscosity and diffusivity respectively after the elimination of n shells, then the 
elimination of the next shell yields the following equations to the first order in perturbation: 

(-iw + v (n) k 2 + 5v (n) k 2 )uf{k) = -^P ijm (k) J dpuf(p)u<(k-p) (9) 
(-iw + H( n )k 2 + 5n( n )k 2 ) ^<(fc) = -ikj J dpuf{p)tl) < {k - p) (10) 



with 



where 



^(n){k) = - {J ^ L2 X ^ * (11) 




5n (n) {k) = ~ Q (12) 



In the above Feynmann diagrams, the solid, wiggly (photon), and curly (gluon) lines represent correlation function 
(uiUj), and Green functions G u , respectively. The filled circle represents (~i/2)Pij m vertex, while the empty 
circle represents — ikj vertex. The RG procedure adopted here is the same as that of Zhou and Vahala [4]. Some of 
the notation used here is close to the that of MHD turbulence calculation of Verma [5, 15]. 

The frequency dependence of the correlation function are taken as: C u (k,uj) = 2C"(fc)5R(G"(fc, uj)) and C*(fc,w) = 
2C^(k)$l(G^(k, u>)). With this assumption, the expressions corresponding to the above Feynmann diagrams will be 

c ,.x 1 [ A dp C u {q) 

Sv [n) {k) = — — — g / —— I S 1 {k,p,q) — - (13) 

(d - \)k 2 7 P+q =k (27r) d v(n){p)p 2 + v(n){q)q 2 

6K(n){k) = ^y P+q . k (^F^ (fc ' ?3 ' 9) «( n )(^ + ^)(,), 2 (14) 
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with 



Si(k,p,q) = kp ((d - 3)z + 2z 3 + (d - l)xy) (15) 
S2(k,p,q) = kp(z + xy) (16) 

The quantities x,y, and z are defined by 

p • q q • k p • k 

gfc p/e 

The effective viscosity and diffusivity after the elimination of (n + 1) shell are 

(ZA «)(„+!) (fc) = (^,K)(„)(fc) + 5(^,K)(„)(fc) (18) 

The spectrum C u (k) can be written in terms of one-dimensional energy spectrum E u {k) as 

where is the surface area of d dimensional spheres. It is known that E u (k) follows Kolmogorov's spectrum, i.e., 

E u {k) = K u (Il u ) 2 / 3 k- 5 / 3 (20) 

where II is the kinetic-energy flux, and K u is Kolmogorov's constant for fluid turbulence. Using the dimensional 
arguments we find that and K(„) have the following forms: 

(u,K) {n) (k n k') = (K u y/\u u y/ 3 k- 4 / 3 (u,Kr {n) (k') (21) 

with k = k n+ \k'{k' < 1). The large-n limit of the v? n \(k') and «L\(fc') are expected to be universal functions in the 
RG sense. 

We solve for v? n \(k') and n* n ^{k') iteratively using Eqs. (13, 14, 18). We take h = 0.7, and start with constant uf Q -. 
and K* Q y We iterate the process till ^* n+1 ^{k') ~ u * n )^') an< ^ K *n+i)(^') ~ K (n)(k'), that is, till they converge. We 
find that the iteration process converges; the limiting value v* and k* are shown in Table I. 

We can draw many interesting conclusions from the above results. Since the scalar does not appear in the equation 
for u, v* computed here is the same as that obtained for fluid turbulence. In Table I we have listed the renormalized 
diffusivity n* and the turbulent Prandtl number Pr tur b- For d = 3, k* = 0.85 and Pr tur b — v* j ' k* — 0.42. The above 
quantities vary a bit with the variation of h, but they are roughly in the same range. The error in our estimate of the 
parameters is of the order of 0.1. Our results are in the same range as those obtained by Zhou and Vahala [4]. 

We have also carried out the above analysis for higher space dimensions. The calculated n* and Prturb are listed in 
Table I. For large d, v* w n* cx d -1 / 2 . The d dependence is in the agreement with the finding of Fournier and Frisch 
for fluid turbulence [16]. The above result also implies that Pr tur b ~ 1 for large d. 

In two-dimensions the scalars are not constrained to double energy-enstrophy conservation like velocity field. The 
RG analysis for two-dimensional scalar turbulence is beyond the scope of this paper. 



III. CALCULATION OF CASCADE RATES 



In this section we compute cascade rates of u and tp, and Bachelor's constant. To this end we use the flux formulas 
and the renormalized parameters computed in the previous section. The time evolution of correlation functions C u 
and (defined by Eqs. [7, 8]) are given by [6, 17, 18, 19] 



where 



(22) 



! +2 ' fc2 ) c " (k ' M) = (a - Dm + *) I +P+q =o (^ [ ^ (k,|p|q) + ^ (k,|q|p)] 

* + **>)c*(M,t) = (27r) ^ k + k0 J k , +p+r _ jgp^WM + ^(k'|q|p)] (23) 



S«"(k'|p|q) = -S([k'.u(q)][u(k').u(p)]) (24) 
S^(k'|p|q) = -3([k'.u(q)][^(fc')V'(p)]) (25) 
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Here S stands for the imaginary part of the argument. Note that Eqs. (22, 23) have been discussed in the earlier 
literature, e.g., Lesieur [6] and Stanisic [17]. However, reinterpretation of the terms S'(klplq) by Dar et al. [19] as 
energy transfer from mode p (the second argument of S) to k (the first argument of S) with mode q (the third 
argument of S) as a mediator makes the formalism more transparent and simple. Also, some quantities which were 
impossible to calculate in earlier formalism could be computed now [19]. This interpretation of Dar et al. is consistent 
with the earlier formalism. 

The energy fluxes U u and from a wavenumber sphere of radius k n is [19] 

Note that there is no cross-transfer between u and tp energy. It is also important to note that both C u and are 
conserved in every triad interaction, i.e., 

S* tlu (k'|p|q) + 5 l,tt (k'|q|p) + S' utt (p|k'|q) + 5 utl (p|q|k') + 5' tI "(q|k'|p) + S ,uu (q|p|k') = (28) 
S'^(k>|q) + 5'^(k , |q|p) + S'^(p|k'|q) + S , ^(p|q|k') + S' ,/ ' ,/ '(q|k , |p) + 5'^(q|p|k') = (29) 

These are the statements of "detailed conservation of energy" in triad interaction (when v = k = 0) [6]. 

The energy fluxes can be calculated using Eqs. (26, 27) by taking ensemble averages of S uu and S^. It is easy to 
check that (S uu ) = (S^) = to the zeroth order, but are nonzero to the first order. The field-theoretic calculation 
performed here is very similar to Verma's MHD flux calculation [18]. Please refer to Verma's paper [18] for further 
details. The Feynmann diagrams for the first order of (S) are 



(S UU (k'\p\q)) = 



(S^{k'\p\q)) = 




(26) 
(27) 



In the above Feynmann diagrams, the solid, dashed, wiggly (photon), and curly (gluon) lines represent (uiUj), 
{tpip), G u , and respectively. In all the diagrams, the left vertex denotes k i: while the filled circle and the empty 
circles of right vertex represent (—i/2)Pij m and — ikj respectively. Algebraically, 

(S uu (k\p\q)) = [ dt'[T 1 (k 7 p,q)G u (k,t-t')G m (p,t,t')G m (q,t 1 t') 



+T 2 (k,p, q)G u (p, t - t')C u (k, t, t')C u (q, t, t') 

+T 3 (k,p, q)G u (q, t - t')C u (k, t, t')C u (p, t, t')] (32) 
(S^(k\p\q)) = f dt'[T i (k,p,q)G^(k,t-t')C^(p,t,t')C u (q,t,t') 

J — OO 

+T 5 (k,p, q)G*{jp, t - t')C^(k, t, t')C u (q, t, t')] (33) 

where Ti(k,p, g)'s are given by 

Ti(k,p,q) = -kp((d-3)z+(d-2)xy + 2z 3 + 2xyz 2 + x 2 z) (34) 
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T 2 {k,p,q) = kp ((d - 3)z + (d - 2)xy + 2z 3 + 2xyz 2 + y 2 z) 
Tz{k,p, q) = kq (xz - 2xy 2 z - yz 2 ) 
T 4 (k,p,q) = k 2 (l-y 2 ) 
T$(k,p,q) = -kp(z + xy) 
We assume the relaxation time for C u {k) and C^{k) to be (^(fc)fc 2 ) -1 and («;(fc)£; 2 ) -1 respectively, i.e., 

C u (k,t,t') = cxp (-v(k)k 2 (t-t')) C u (k,t,t) 
C^{k,t,t r ) = exp(- K (/c)fc 2 (i-t'))C" /J (fc,t,i) 
With this assumption, Eqs. (32, 33) reduce to 



(35) 
(36) 
(37) 
(38) 

(39) 
(40) 



n"(*o) 



n*(*o) 



dp 



k>k W J P <k i^Y v(k)k 2 + y{p)p 2 + v{q)q 2 
[Tifap, q)C u (p)C u {q) + T 2 (k,p, q)C u {k)C u (q) + T 3 (k,p, q)C u {k)C u (p)] 
dp 1 



Jk>k (27r) d J p< 



J p<ko (2ir)d K (fc)fc2 + K ( p )p2 + „( g ) g 2 
[T 4 (k,p,q)C^(p)C m (q)+T 5 (k,p,q)C^(k)C m (q)} 
For C u {k) we substitute Eqs. (19, 20), while for we substitute [6] 

C*(jfe) = &l k -^E u (k), 

E^(k) = i^ir^rp)- 1 / 3 *;- 5 / 3 

where is called the Batchelor's constant. The renormalized viscosity and diffusivity in the inertial range are 



(K u ) 



1/2(^1/3^-4/3^ 



K (jfe) = (^) 1 / 2 (n u ) 1 / 3 fc- 4 / 3 K *. 

The substitution of the above quantities, and the change of variables 

fc fco fco 

fc = — ;p = — v. q = — w 

u u u 

yield the following nondimcnsional version of the flux equations [16]: 

, r aq. - r 1 r 1 + v 
(K u ) 3/2 



i— / dv\n(l/v) / rfw(i;u;) d - 2 (sina)' i ~ 3 f tl (w, 



1 



(d-1) 2 



w) 



(d- 



l+u 



(iw(™) <i - 2 (sina) <i - 3 F' , '(i ) , W ) 



1— V 



where a is angle between vectors p and q, and the integrals F u '^(v, w) are 



pu 



1 



u*(l 



,2/3 



2/3) 



[ti(v,w)(vw) d 3 +t 2 (v,w)w d 3 +t 3 (v,w)v d 3] 



K*(l + W 2 / 3 ) + J/*UI 2 



■j^[U(v,w)(vw) d 3 +t 5 (v,w)w 



(41) 

(42) 

(43) 
(44) 

(45) 
(46) 

(47) 

(48) 
(49) 

(50) 
(51) 



Here ti(v, w) — Ti(k,kv, kw)/k 2 . 

The terms in the square brackets of Eq. (48, 49) (denoted by I u ^) involve integrals. We compute them using 
gaussian quadrature. The integrals converge for all dimensions d > 2. Once the integrals are known, Kolmogorov's 
and Batchelor's constants (K u and respectively) can be computed. The computed values are given in Table 1. 

In our calculation Batchelor's constant in three dimension is 1.25. Due to uncertainties in the value of v* and 
k* , the error in the constant could be of the order of 0.1. Earlier, Kraichnan had estimated the constant to be 0.2. 
Yakhot and Orszag [2] obtained = 1.16 by their e-based renormalization group analysis. Canuto and Dubovikov 
[11] and Canuto et al. [12] estimated = (5/3) * 0.72 = 1.2 using their RG calculation. Lin et al. [7] find the 
constant to be close to 0.3. Our result is in very good agreement with the theoretical predictions of Yakhot and 
Orszag [2] and Canuto et al. [12], as well as to the experimental values (« 1.2 — 1.4, sec Monin and Yaglom [20]). 

It is also interesting to note that both K u ^ are proportional to d 1 / 3 , consistent with the predictions of Fournier 
and Frisch [16] for fluid turbulence. This result implies that the cascade rated IP''* will decrease with dimensions as 

d-V2. 
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TABLE I: The computed values of renormalized viscosity v* , diffusivity n* , turbulent Prandtl number Pr tur b, Kolmogorov's 
constant K u and Batchelor's constants for various space dimensions d. 



d t£ Pr turb FT 

3 0.36 0.85 0.42 1.53 1.25 

4 0.42 0.69 0.61 1.60 1.39 
7 0.38 0.48 0.80 1.76 1.65 
10 0.34 0.39 0.87 1.94 1.83 
25 0.22 0.24 0.94 2.43 2.44 
50 0.16 0.16 1.0 3.1 3.0 
100 0.093 0.095 0.98 3.4 3.4 



IV. CONCLUSIONS 

In this paper we employed field-theoretic techniques to calculate the cascade rates of scalar turbulence. Our 
calculation is to first order. From this formalism we also calculate Batchelor's constant. In three dimensions, we 
find Batchelor's constant to be 1.25, which is in very good agreement with the theoretical predictions of Yakhot and 
Orszag [2] and Canuto et al. [12], and the experimental values. In higher space dimensions the constant varies as 

Our calculation of cascade rate requires the renormalized viscosity and diffusivity. We have extended the RG 
calculations of Zhou and Vahala [4] for higher dimensions. Our calculations show that for higher dimensions, the 
renormalized viscosity and diffusivity vary with dimensions as dr 1 !" 1 , and the turbulent Prandtl number approaches 
unity. 
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